Spin Hall Conductivity on the Anisotropic Triangular Lattice 
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We present a detailed study of the spin Hall conductivity on a two-dimensional triangular lattice 
in the presence of Rashba spin-orbit coupling. In particular, we focus part of our attention on the 
effect of the anisotropy of the nearest neighbor hopping amplitude. It is found that the presence 
of anisotropy has drastic effects on the spin Hall conductivity, especially in the hole doped regime 
where a significant increase or/and reversed sign of the spin Hall conductivity has been obtained. 
We also provide a systematic analysis of the numerical results in terms of Berry phases. The changes 
of signs observed at particular density of carriers appear to be a consequence of both Fermi surface 
topology and change of sign of electron velocity. In addition, in contrast to the two-dimensional 
square lattice, it is shown that the tight binding spin-orbit Hamiltonian should be derived carefully 
from the continuous model on the triangular lattice. 

PACS numbers: 72.25.Dc,73.23.-b,71.20.-b,71.70.Ej 



I. INTRODUCTION 

The incorporation of electron spin degree of freedom 
and its associated magnetic moment into electronic de- 
vices has initiated a rapidly growing field of spin elec- 
tronics, called spintronics. One of the first idea of sphi- 
tronic device is attributed to Datta and Das in 1990 1 . 
The authors have proposed a spin field-effect transistor 
with the innovation of manipulating a pure spin current. 
Such a manipulation is possible due to the rclativistic 
effect of the coupling between spin and orbital degrees 
of freedom which can be manipulated e.g. by means 
of a gate voltage. This spin-orbit coupling (SOC) can 
emerge in various ways in semiconductors. One source 
is the lack of inversion symmetry of the confining po- 
tential, which defines the 2D electron gas, generated by 
the hcterostructure of the semiconductor. This Bychkov- 
Rashba 2,3 type SOC has the advantage that it can be 
directly tuned by an applied gate voltage. Another way 
to introduce SOC in the sample is by choosing a noncen- 
trosymmetric material where the lack of inversion cen- 
ter creates a SOC called Dressclhaus 4 SOC. Developing 
further the correspondence between (charge) electronics 
and spintronics with the aim of creating pure spin cur- 
rents, D'yakonov and Perel 5 proposed already in 1971 
the spin Hall effect (SHE). This mechanism requires spin 
dependent impurity scattering such as skew scattering 
or side jump mechanism 6 . In today's terminology, it is 
known as extrinsic SHE. It was experimentally confirmed 
in 2004/2005 by angle-resolved optical detection of spin 
polarization at the edges of a two-dimensional layer. 7,8 

In this paper, we will focus on the SHE which arises 



even in the absence of impurities, the so called intrinsic 
SHE, which is due to SOC. The latter is assumed in 
this paper to be of Bychkov-Rashba type. It lifts, 
without breaking time reversal symmetry, the spin 
degeneracy of the eigenstates. The theory of the SHE 
has been developed in the last ten years, as reviewed in 
Refs. 6,9, but only recently experimental evidence has 
been reported for electrical manipulation and detection 
of intrinsic SHE in ballistic HgTe/HgCdTc quantum 
wells. 10 

Starting at low electron fillings in the conduction band 
where both Rashba bands are occupied, one finds the 
"universal" spin Hall conductivity (SHC) ctsh = e/8ir, 
which is independent of the SOC strength for sufficiently 
small couplings. 11,12 Going beyond the long wavelength 
approximation, the underlying lattice geometry becomes 
important. On a square lattice using a single-orbital 
model (spin-split) the SHC ctsh shows electron-hole 
symmetry: asu is an odd function of the Fermi energy 
Ep which thus vanishes at half-filling. 11 ' 12 However 
both, shape and sign of SHC are lattice structure 
dependent. This can, for instance, nicely be illustrated 
with the appearance of several plateaus as seen in the 
SHC calculations on honeycomb and kagome lattice by 
Liu et ai, Ref. 13,14. Both of these lattices contain 
several atoms per unit cell respectively two and three. 
The kagome lattice can be seen as a triangular lattice 
with 3 atoms per unit cell. The question which arises 
is whether similar differences as seen between square 
and kagome SHC are expected already in the case of 
the simple triangular lattice. In contrast to the square 
lattice, the triangular lattice brings along the property 
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of geometrical frustration which is at the origin of 
exotic properties. In addition, electron-hole symmetry 
is broken on the triangular lattice. Among exciting 
properties, the triangular based lattices are known to 
favor in some cases ferromagnetic groundstates when 
the electron-electron interaction is switched on. 15-17 
This is also known as flat band ferromagnetism. 18-21 
In these systems, the intrinsic frustrating nature of the 
lattice is a crucial ingredient in the stabilization of 
ferromagnetism. The importance of geometrical frustra- 
tion has also been shown to play an important role in 
supcrfluid-Mott insulator quantum phase transitions in 
ultracold quantum gases. 22 It has been suggested that 
the supersolid phase reported results from the com- 
petition between Mott localization and frustration. 23 
Among other interesting features found on triangular 
lattices one can also mention superconductivity in C0O2 
layered compounds 24 or various anisotropics in organic 
biscthylcncdithio conductors 25 . 




Figure 1: (Color online) Two dimensional triangular lattice 
in the x-y plane, t and t are the hopping amplitudes, a is the 
lattice constant. The directions are numbered. 



Back to the SHE, it is interesting to mention that there 
are already experimental realizations of surfaces with tri- 
angular geometry. Indeed, giant Rashba SOC has been 
recently reported on surfaces of Au(lll), Bi(lll) and 
Bi/Ag(lll). 26 In this manuscript we present a detailed 
analysis of the SHC on the anisotropic triangular lattice 
as a function of the carrier concentration. In addition, 
the obtained numerical results will be discussed in the 
context of Berry phases 27 . One of the aims of including 
anisotropy in the hopping amplitudes is to allow us to 
tune in a continuous way the position of the Van Hove 
singularity and in the same time change the topology 
of the Fermi surface. Experimentally, the anisotropy in 
the hopping integrals could be achieved by introducing 
strains or growing samples on surfaces which have small 
differences in the lattice parameters. 



As it will become clear in the following, it is convenient 
to introduce the lattice unitary vectors, 



The Hamiltonian can be recast in the following form, 
H =5^((e x ■ k) 2 + (e 4 ■ k) 2 + (e 5 ■ k) 2 ) - m* a\ 



= l^(kl + + fcjj) - y (fc 4 - k + 2k x )(Ty 
+ -^=(fc 4 + h)a x , 

where the canonical moment is given by 



(6) 
(7) 



II. TIGHT BINDING HAMILTONIAN ON 
TRIANGULAR LATTICE WITH SOC 

We consider a 2D electron gas on the triangular lattice 
as sketched in Fig. 1. 

Including the Rashba SOC with coupling strength a 2 
the Hamiltonian reads 



H 



2m 
1 



+ a 2 (k y a x - k x Oy) 



-(k + eAso) -m*a 2 , 



where, 
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(2) 

(3) 
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and <Tj, i = x,y are the Pauli matrices, S is the spin- 1/2 
operator and m* the electron effective mass. 



k =k + a 2 m* 



(8) 



The discretization of the Hamiltonian leads to the Tight 
Binding Hamiltonian (TBH) whose nearest neighbor hop- 
pings are listed in Tab. II. This Hamiltonian is identical 
to that one would obtain directly using the well known 
tight-binding expression of the SO part and that reads, 



H (R) = 2 «2 

3 2a ^ 



(ij)aa' 



(T X 



(9) 



where ey is the unitary vector pointing from site i and j. 
Note the presence of the 2/3 coefficient in both Eq. (7) 
and in Eq. (9). We could absorb this coefficient in the 
hopping definition but we have chosen to keep them in 
order to facilitate the direct comparison to existing calcu- 
lations. The comparison between our TBH and that de- 
rived in Ref. 28 reveals differences in the hopping ampli- 
tudes. The difference in the derived TBH results from the 
discretization procedure. By using the relation between 



3 



{hj) 
(0,1) 
(0,2) 
(0,3) 
(0,4) 
(0,5) 
(0,6) 



3 1 ^/%° x % 3 a y 



10 



in 



R . tR 



"3*0 -« 75 <7 a! + «- 



3 a V 



3^0 
2 * 



7T 



1 3 a v 



Table I: Matrix elements in the tight-binding model for the 
kinetic part and the Rashba part H^ R ' , where t R — 

02/ (2a), to — l/(2m*a 2 ) and a denotes the lattice constant. 



k x , k y and ki, k§ one can generate an infinite number 
of rigorously equivalent Hamiltonians in the continuous 
picture with different coefficients for a x and a v . How- 
ever, the discretization will lead to an infinite number of 
inequivalent TBHs. Thus this crucial step has to be done 
carefully. Note that this problem will not occur on the 
square lattice because the connecting vectors e x and e y 
are orthogonal to each other and the discretization can be 
done unambiguously. Then, how to proceed in order to 
derive the correct TBH? There are two procedures. The 
most straightforward is to use directly Eq. (9) as often 
done in the literature. The second procedure consist in 
re-expressing the continuum Hamiltonian in term of the 
canonical momentum as done in Eq. (7) and then only 
discretize it by properly defining the real space deriva- 
tive. 



Spectrum and Density of States 
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Figure 2: (Color online) Degeneracy points in the energy 
spectrum of H in fc-space (for the 1st Brillouin zone shown): 
M = Pi = {0,-2tt/^3} T , Pi = {2tt/3, -2tt/^3} t , P s = 
{7r,-7r/^} T , K = P 4 ^ {471-/3, 0} T , P 5 = {7r,7r/v^} T . The 
blue dots indicate the saddle points in case of vanishing SOC. 



The corresponding eigenvectors are 



|A ± (k)) 



1 / ±e i0(k) 



V2 



1 



(12) 



with </>(k) given by 



e ^(k) = ^=*i(k) ( cos 
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+ sin(A; x ) + i\/3cos ( sm f 



(13) 



Let us now proceed with the diagonalization of the 
anisotropic TBH. We define the hopping amplitude in 
the x and — 4 direction as t = (2/3)frj and that in the 
— 5 direction as t = (2/3)to- From now on, the lattice 
constant a will be set to one, a = 1. 
The spectrum of the TBH is given by 



E±(k) = - - (2t (cos(k x ) + cos (fc 4 )) 
+2t cos (k 5 ) ± v^^-F^k) 



where -Fi(k) is defined as 

Fi(k) = (3 + cos(fc x ) - cos(2fc x ) 

— (1 + 2 cos(fc x )) cos(v / 3fc y ) 

VSky \ 
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Before proceeding with the calculation of the SHC, let us 
first discuss some features of the energy spectrum which, 
as will be seen, will appear to be important and help- 
ful to understand the characteristics of this dynamical 
quantity. As it is well known, the SHC can be directly 
related to the fc-space curvature, 29 this will be presented 
in detail in what follows, Sec. III. In this context, the 
singular values of the Berry connection at degeneracy 
points of the energy spectrum are of crucial importance. 
To prepare the insight into the topological nature of this 
problem, we first calculate the degeneracies by solving 
directly Fi(k) = 0. The result is plotted in Fig. 2. 

In the absence of SOC (t R =0) and in the isotropic case 
(to = to), the density of states (DOS) exhibits a Van Hove 
singularity at Fermi energy E = 4/3to- It corresponds to 
a straight line which goes through P3 and P5. For finite 
values of SOC strength, in the isotropic case, this peak 
is now split as clearly seen in Fig. 4. In addition, one 
has two additional peaks at the boundary of the energy 
spectrum. 30 The energy dispersion is plotted as a func- 



(11) tion of momentum in Fig. 3(d). As seen, the degeneracy 
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Figure 3: (Color online) The energy dispersion, Eq. 10, is plot- 
ted for both modes E± (red curve: E-). The fe-points P; are 
according to Fig. 2. The two red horizontal lines indicate Esi 
and Es2- The green-dashed one indicates Edi and the orange- 
dashed Edi- We have used a large SOC coupling, t R — lfo, 
for sake of clarity. (See text) 



points correspond to the points Pi (i = 1, 2, . . . , 5) and T. 
One finds E{P 1 ) = E(P 3 ) = E{P 5 ) and E(P 2 ) = E(P A ). 
In Fig. 4 these particular energies arc indicated by dashed 
lines. When we switch on the anisotropy the singularities 
in the DOS change significantly as expected. In Fig. 3 one 
now sees that one degeneracy is lifted E(P{) = E(P 3 ) ^ 
E(P 5 ), and we still have E{P 2 ) = E(P 4 ). Thus in the 
presence of anisotropy we have four relevant energies that 
we denote, 





E(T) 


= Edi 


= -(4/3)(2 + toAo) 




E(P 5 ) 


= Edi 


= (4/3)(2-? /to), 


E{Pi) 


= E(P 3 ) 


= Egi 


= (4/3)t /to, 


E(P2) 


= E(P 4 ) 


= Es2 


= (4/3) + (2/3)t /io, 



(14) 



with the energy expressed in values of io- Note that, 
two different sets of labels (di,^) and (Si, S2) are intro- 
duced. The meaning of this separation will become clear 
in what follows. 



III. SPIN HALL CONDUCTIVITY 

Numerical Results of SHC 

Having analyzed the energy spectrum, we now pro- 
ceed further with the calculation of dynamical properties, 
more precisely the spin Hall conductivity ctsh- The SHC 
is calculated within linear response theory using Kubo 




Figure 4: (Color online) DOS p(E) for a system of size L — 
300 (cutoff in the Lorenzian was chosen to be 77 = 0.03) and 
t R = 0.3to, for different values of to/to = 0, 0.2, . . . , 1. Dashed 
lines indicate the degeneracy points energy, see Fig. 3. 



formula 
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asn(E F ) = 2 



V 



E 

E m <E F <E, 



S((V> m | Jl |V> n ) (lp n \ V y \lp m )) 

(E n - E m y + Ty 2 



(15) 

where V is the total number of lattice sites, Ei the 
cigenencrgies, \ipi) the corresponding cigenstates, Ep the 
Fermi energy and rj a positive infinitesimal. The spin 
current operator is 



(16) 



and the velocity is defined with the usual expression v = 
i [H, r] . To calculate the SHC wc first compute the matrix 
element density function j(x, y) defined by 

j(x, y)=J2 M mnS(x - E m )S(y - E n ), (17) 



(18) 



where the matrix element M mn is 

M mn = $S[{lp m \ Jl \lp n ) (lpn\v y \lp m )}- 

Note that the function j(x, y) is a particularly useful and 
convenient quantity for calculations carried out by Kernel 
Polynomial Method 34 in the presence of disorder. 12,35 

We can thus rewrite asa(Ep) in terms of j(x, y), lead- 
ing to the following expression 12 : 

osn(E F ) = 

(19) 



e 
V 



(y — x) 2 + rf 



'_E mi „ ^B min 

where Je f is the Fermi-Dirac distribution. 

In the present case the only non-zero-matrix elements 
are 

M±(k) = - M T (k) (20) 

= 3{(A+(k)| Jl |A"(k)) (A-(k)| vy |A+(k))}, 

(21) 
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(a) 



Figure 5: (Color online) Matrix element density function 
j(x,y) for the isotropic case (to = to) [(b) is a 2D plot of 
(a)]. The system size is of 70 x 70 sites. The SOC strength is 
t R — O.Mo- The dashed lines in (b) mark the position of Esi- 
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Figure 6: (Color online) (a) SHC for different system sizes 
L = 50, . . . , 300 for the isotropic case. The SOC strength 
is t R — 0.3to- (b) Comparison between the ctsh calculated 
within our model (continuous line) and that of Ref. 28 (filled 
circles). The system contains 70 x 70 sites. The vertical 
dashed line corresponds to .Esi. 



where 

M±(k) = 4V2i*sin(fcx) + 2 cos (^j cos Mj^ J 
x (2i sin(k x ) + t sin (fc 4 ) - t sin (& 5 )) /(3Fi(k)). 

(22) 

The function j(x,y) is plotted in Fig. 5 for the isotropic 
case. Notice the symmetry due to Eq. 20. One can 
clearly see a sign change for j values at (x, y > x) resp. 
(x, y < x). This gives a clear indication for a sign-change 
in osh at a particular Fermi level. In fact, this energy co- 
incides with Esi- The resulting SHC is plotted in Fig. 6. 
At low electron density we observe a sharp increase to 
csh = e /(87r) as expected since the lattice structure is 
irrelevant. As we increase the Fermi energy further, cgH 
remains almost flat up to E/to ps 1 (extended plateau). 
Then at E = Esi the sign of SHC changes in agreement 
with the discussion on j(x,y). More remarkable is the 
change of amplitude of the SHC in the region of negative 
conductivity. Indeed, it increases by 60% with respect to 



the case of low electron density. 

Interestingly, the calculations indicate that the finite 
size effects are almost negligible up to Esi, but not in 
the hole doped regime. As mentioned before, one has to 
be careful in the derivation of TBH. We now compare 
the calculated SHC with that obtained using the TBH of 
Ref. 28. The result is plotted in Fig. 6 (b). As seen, the 
results are relatively comparable up to Ep = —2to and 
drastically deviate in the vicinity of Esi ■ The conductiv- 
ity calculated within the TBH of Ref. 28 exhibits a clear 
maximum in the vicinity of Esi, but no change of sign. 
We now analyze the effect of anisotropy in the hopping 
amplitude on the SHC. The results are depicted in Fig. 7. 
As we discussed it previously a change of sign occurs at 
Esi- As we reduce to we observe a shift of the cross- 
ing energy towards Ep = 0. One also notices drastic 
changes in the shape of <tsh beyond Esi- In particu- 
lar, we observe two minima and a significant reduction 
of the amplitude of SHC. This is especially pronounced 
for the case to = 0-7t - As the anisotropy is further de- 
creased, the second minimum disappears and a second 
change of sign of the SHC is found. This occurs at the 
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Figure 7: (Color online) SHC for t R = 0.3f and to /to = 
0, 0.1, . . . , 1. The local maximum at low energies is shifting 
with Edi- The sign-change at low Fermi energy is found at 
Esi, the other one (in the case of to /to < 2/3) is found at 
Es2- The vertical dashed line corresponds to -Esi (to = to) = 
E S2 (to = 0). 



value of to = (2/3)fo. Thus, below this value one has 
two changes of sign, one at Esi and the second at E$2- 
One natural question which arises is why beyond this 
value SHC changes its sign a single time? This will be 
answered and understood when discussing the SHC in 
the framework of Berry phases. Remark that the second 
change of sign could not be anticipated from the DOS. 
Indeed, in contrast to the peak seen at Esi no peak, no 
singularity or pronounced feature is visible at Es%, see 
for instance the case to = (l/3)frj ln Fig. 4. The limiting 
case to = requires also some additional attention. In 
the absence of SOC, the system is topologically equiva- 
lent to the square lattice. As seen in Fig. 7 the results for 
<7sh are similar to that of the square lattice for Ep < 
(see e.g. Rcf. 12,36,37). However, for Ep > (below 
Ep/to < 1), the amplitude is reduced to half of that 
of the square lattice. The spin flip term in the (2 — 5) 
direction at finite SOC is at the origin of these crucial 
differences in the hole doped regime. It is surprising that 
the drastic changes take place only for Ep > 0, we will 
shed light on it in what follows. 



Geometric Interpretation of SHC 

In order to get a better understanding of the previously 
discussed features, we propose to analyze the results in 
the context of the Berry phase. As could be guessed from 
comparing the spectra shown in Fig. 3 and the SHC plot- 
ted in Fig. 7, the topology of the spectrum plays an im- 
portant role in understanding the transport properties of 
electrons and their spins. We remind that the transport 
of a Bloch state |A(C)) in the parameter space C which 
has a non- vanishing curvature leads to picking up a geo- 
metrical phase. In the case of adiabatic transport, valid 



6 
4 
2 

*r o 

-2 
-4 
-6 



-y//o) 



-6-4-2 2 4 6 

kx 

Figure 8: Berry connection A ± (k) shown as a stream plot 
in k-space. Divergences of the Berry connection A are found 
at the degeneracy points of the energy spectrum indicated 
by small circles. Clockwise/counterclockwise winding of A 
around this degeneracies is indicated by yellow/black color of 
this circles. 



in the linear response regime, it is called Berry phase. 27 
In momentum space this geometrical phase is at the ori- 
gin of the quantum, 38-40 the anomalous 41 and spin 42 ' 43 
Hall effects. The latter was also studied in the context 
of Fermi-surface topologies on honeycomb and kagome 
lattice by Liu et al. in Ref. 13,14. 

In our case the Berry phases are explicitly given by the 
integral 
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Jc± 



(23) 



along a loop C = {k x ,k y } in momentum space with the 
condition E ± (k x , k y ) = Ep, the Berry connection ^l ± (k) 
is 



A ± (k)= (A ± (k) 







idk 



A ± (k) 



(24) 



In our case we have found that -4 + (k) = _4~(k). In- 
deed, the eigenvectors, Eq. 12, do not depend on the hop- 
ping anisotropy and therefore the Berry connection is in- 
dependent of to /to- The Berry connection is shown as 
a stream plot in Fig. 8. Examining it in the vicinity of 
the degeneracy points of the energy spectrum, T and Pi 
(i = 1, . . . , 5), we find different winding directions for A. 
It diverges exactly at these points. In order to relate our 
findings to the Berry connection we reformulate the SHC, 
Eq. (15), in terms of it, 14,44 



ctsh(^f) = I 



[A m Ck) 



(k)] 



(25) 
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Figure 9: (Color online) Fermi surface for both spin-split 
bands (blue: E+, according to Fig. 3) for (a) Ef = Esi — e 
and (b) Ef = Esi + e for to = to and t R = O.lio. The arrows 
indicate v ~ v' ' . 



where v^°^(k) is the velocity in absence of SOC. 
Let us show that the SHC sign is connected to that of 
[>4 m (k) x v(°)(k)] z . For that purpose we first consider 
the isotropic case. Notice that below Esi, as the term 
(E + (k) — £L(k)) is always positive, the sign change of 
SHC is associated to that of [A m (k) x V (°)(k)] z . For 
Ep < Esi the winding direction is always clockwise for 
the Berry connection and the velocity is always pointing 
outside (see Fig. 9 (a)), resulting then in a positive SHC. 
Beyond Esi the Fermi surface topology is completely dif- 
ferent, it consists now of closed loops around the P2 and 
P4 as illustrated in Fig. 9 (b). Furthermore, the direc- 
tion of the velocity is now directed inside (towards Pi 
and P4). However, the winding directions (see Fig. 8) of 
the Berry connection is unchanged, clockwise as around 
the r point. This results in a sign change of the spin 
Hall conductivity. We now proceed further and discuss 
the presence of anisotropy. Since as mentioned before, A 
is independent of to thus the changes in SHC are associ- 
ated with the Fermi surface topology only. Let us start 
with the extreme case of a fully anisotropic system, e.g. 



to = 0. As seen before, the SHC exhibits two changes of 
sign, one at Esi and the other at E$2- Thus, we focus 
our attention on these two Fermi energies. The situa- 
tion at Esi, Fig. 10 (a) and (b), is very similar to that of 
the isotropic case, except the fact that the Fermi surface 
at Ep = Esi + e is now also surrounding the point P5 
for which the winding direction of Berry connection is 
now negative (anti-clockwise). As the velocity direction 
is still pointing inwards, the negative winding around P5 
can lead to a sign change in [„4 m (k) x v^°'(k)] z , thus, to 
the appearance of terms in Eq. 25 which reduce the ab- 
solute value of <7sh- This explains why osh varies from 
~ — 1.6e/(87r) in the isotropic case to ~ — 0.5e/(87r) for 
to = . Concerning the second sign change, we find at 
Ep = Es2 ± e the same velocity direction. However, the 
Fermi surface at Ep = Es2+£ is, in contrast to Fig. 10 (b) 
and (c), surrounding only P5 which as mentioned before 
has the opposite winding direction with respect to that 
of P2 and P4. This yields the topological explanation of 
the change of sign at Es2- 

IV. CONCLUSION 

In this paper we have performed a detailed theoreti- 
cal study of the spin Hall conductivity on the anisotropic 
triangular lattice in the presence of Rasbha spin-orbit 
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coupling. Note that experimentally anisotropy could be 
e.g. tuned by taking advantage of lattice mismatches 
in semiconductor heterostructures. To conclude, it has 
been shown that on such a lattice the tight binding model 
should be carefully derived from the continuous Hamil- 
tonian. In addition, our calculations have revealed that 
anisotropy in the hopping term has drastic effects on the 
SHC especially in the hole doped regime. Furthermore, 
depending on the anisotropy strength, one finds either 
a single or two changes of sign in the SHC at particu- 
lar carrier densities. The origin of the observed changes 
in the SHC has been interpreted and understood geo- 
metrically in terms of Berry phases. The high symme- 
try points which have different winding directions of the 



Berry connection control the rapid variation of the SHC. 
More specifically, the changes of sign and amplitude of 
the SHC have been associated to the topology of the 
Fermi surface, the number of high symmetry points en- 
closed in it and to the sign of the velocity in the absence 
of SO coupling. 
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